Ground-State Energy and Condensate Density of a Dilute Bose Gas Revisited 
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The ground-state energy per particle E/N and condensate density no of a dilute Bose gas are 
studied with a self-consistent perturbation expansion satisfying the Hugenholtz-Pines theorem and 
conservation laws simultaneously. A class of Feynman diagrams for the self-energy, which has 
escaped consideration so far, is shown to add an extra constant d p ~ 0(1) to the expressions 
reported by Lee, Huang, and Yang [Phys. Rev. 106 (1957) 1135] as E/N = (2Tvh 2 an/m)[l + 



(128/15y'¥ + 16ci p /5)%/a 3 n] and no/n = 1 — + Ci p )Va 3 n, where a, n, and m are are the 

s-wave scattering length, particle density, and particle mass, respectively. We present a couple of 
estimates for Ci p ; the third-order perturbation expansion yields Ci p = 0.412. 
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Lee et alX studied the ground-state properties of a 
single-component homogeneous dilute Bose gas based 
essentially on Bogoliubov theory^ at zero temperature. 
They obtained widely accepted expressions for the en- 
ergy per particle and condensate density as-^ 1 
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where a, n, and m are are the s-wave scattering length, 
particle density, and particle mass, respectively. On the 
other hand, a standard finite-temperature extension of 
Bogoliubov theory is known to yield an unphysical en- 
ergy gap in the single-particle spectrum^ in contradic- 
tion to the Hugenholtz-Pines theorem.— This issue was 
discussed extensively by Hohenberg and Martin* in con- 
nection with conservation laws. 

Recently, a self-consistent perturbation expansion has 
been formulated for condensed Bose systems so that 
the Hugenholtz-Pines theorem and conservation laws are 
satisfied simultaneously^^ A notable prediction of this 
theorjjii is that there is a new class of Feynman diagrams 
for the self-energy i.e., those that may be classified as 
improper or reducible following the conventional termi- 
nology and should certainly be excluded from its defini- 
tion in the normal state^^^ The rationale for their ex- 
istence in condensed Bose systems is that they are indis- 
pensable for three exact statements, i.e., the Hugenholtz- 
Pines theorem, conservation laws, and an identity for the 
interaction energy, to be satisfied order by order simul- 
taneously in the self-consistent perturbation expansion. 9 
Thus, the conventional procedure of constructing the self- 
energy in the Dyson-Beliaev equation with the concept of 
properness (or irreducibility) 3 ,4 cannot be justified. Since 
they have been overlooked, it will be worth looking at 
how the improper diagrams affect the standard results 
for the dilute Bose gasi^ 

In the present paper, we will show that the improper 
diagrams change eq. (JTJ into 
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where c lp is a numerical constant of order 1 . We will sub- 
sequently estimate c; p using a couple of approximations. 
One of them, i.e., the perturbation expansion up to the 
third order, will be shown to yield q p = 0.412. 

We consider a homogeneous system of identical Bose 
particles with mass m and spin interacting via the con- 
tact potential U8{r\ — r 2 ). The Hamiltonian is given by 
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where p, e p , \i, and V are the momentum, kinetic en- 
ergy, chemical potential, and volume, and c£ and c p are 
the creation and annihilation operators, respectively. It 
is convenient to set h = 2m = fee = T~ = 1, where &b 
denotes the Boltzmann constant and T c ° is the transition 
temperature of ideal Bose-Einstein condensation. 4 Thus, 
the kinetic energy is expressed simply as e p = p 2 . Ultra- 
violet divergences inherent in the continuum model are 
removed here by introducing a momentum cutoff p c 3> 1 . 
However, our final results will be free from p c , as seen 
below. It is standard in the low-density limit to remove 
U in favor of the s-wave scattering length a. They are 
connected in the conventional units by 
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with 0(x) the step function, which in the present units 
reads l/87ra = 1/U + p c /4ir 2 . We will focus on the limit 
a — > and choose p c so that 1 <C p c <C 7r/2a is satisfied. 
Thus, we can set 
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as an excellent approximation for the present purpose. 

Let us recapitulate some relevant results from ref. 0. 
Green's function for the condensed Bose system can be 
expressed in the Nambu representation as 



Gp = 



—Gp* 



(5a) 
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where p = (p, ie e ) with e e = 2£ttT (I = 0, ±1, ±2, • • • ) 
the Matsubara frequency. The upper elements satisfy 
Gp = Gp» and Fp = F_p, and a barred quantity generally 
denotes Gp = G*_p» . The matrix Green's function obeys 
the Dyson-Beliaev equation 
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which may also be regarded as defining the self-energies 
Yip and A^. In the self-consistent perturbation expan- 
sion, they are obtained from a functional $ [Gp, uq] as 



^'be- ^-^w.- (6a) 



In addition, $ satisfies 
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The substitution of eq. (|6a|) into eq. (|5bj) yields self- 
consistent (i.e., nonlinear) equations for Gp and Fp. It 



also follows from eq. (|6b[) that the stationarity condition 
dfl/SriQ — for the thermodynamic potential f2 is equiv- 
alent to the Hugenholtz-Pines relation 
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These are exact statements. It has been shown that the 
key functional $ can be constructed as a power-series 
expansion in U in such a way that eq. ([7]), conservation 
laws, and an exact relation for the interaction energy are 
fulfilled simultaneously order by order. 

The first-order functional $W is given graphically in 
Fig. [U which can be shown to reproduce eq. JT]) as follows. 
The differentiations of eq. (|6aj) correspond to removing 
a line of Gp and Fp, respectively, from every diagram 
in Fig. [T] in all possible ways. Hence, Y^ and are 
obtained as 



= 2Un, 
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where 
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is the particle density with 0+ an infinitesimal positive 
constantj&i^ and the summation over p denotes 
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Thus, the first-order self-energies Y^ 1 * 1 and A^ 1 ) have no 
p dependence, and eq. ([7]) reduces to 
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Substituting eqs. ([8a)l. (|8b|). and ((8d) into eq. (J5bj> and 
carrying out matrix inversion, we obtain the upper ele- 
ments of Gp in eq. (|5a[) as 
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With eq. ([9]) , we can perform the summations over p in 
eqs. (|8b|) and (|8cjl analytically to obtain 
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where we have set p c — > oo in the second convergent 
integral; this procedure will be adopted throughout. The 
substitution of eq. ((4]) into eq. (|8al) yields 



(11a) 



= Wnan ^1 + -p c a^ . 



For AW, eqs. (0} and (JTOJl indicate that A^ sa 87ran to 
the leading order. Collecting next-order terms perturba- 
tively with p c a <C 1 in mind, we arrive at 
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Let us substitute eq. (fTT]) into eqs. (fHH]) and (|10b|) . We 
then obtain 



FIG. 1: Feynman diagrams for A filled circle denotes 

2(7, a line with an arrow (two arrows) represents G (either 
F or F) in eq. ()5a[) as in the theory of superconductivity, 3,4 
and every missing line in the last three diagrams corresponds 
to no. The number below each diagram indicates its rela- 
tive weight, which should be multiplied by 1/4 to obtain the 
absolute weight. 
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and eq. (|lb|). respectively, which are correct up to the 
next-to-the- leading order and also free from p c . Finally 
integrating the thermodynamic relation /i = dE/dN over 
the particle number N , we arrive at eq. (|la[) in units of 
h = 2m = 1. Thus, the Lee-Huang- Yang results of eq. (H} 
have been reproduced in our mean-field approximation. 



FIG. 3: Feynman diagrams for <3>' 3 ' drawn without arrows. 



FIG. 2: Feynman diagrams for $' 2 '. A filled circle denotes 
2(7, a line with an arrow (two arrows) represents G (either F 
or F) in eq. (|5a|l . and every missing line in the second and 
third rows corresponds to no. The number below each dia- 
gram indicates its relative weight, which should be multiplied 
by —1/8 to obtain the absolute weight. 



Now, we proceed to look into higher-order terms. The 
second-order $( 2 ) can also be constructed uniquely so as 
to satisfy eq. ([?]) . conservation laws, and an exact relation 
for the interaction energy^ The results are shown graph- 
ically in Fig. [3J For the dilute Bose gas at T = 0, in gen- 
eral, the relative importance of each diagram decreases as 
the number of Green's function lines in it increases, as our 
previous analysis on eq. © indicates. Hence, the domi- 
nant contribution of the second order originates from the 
third-row diagrams in Fig. [2] which also bring a novel 
structure to the self-energies. Indeed, the differentiations 
of eq. (j5a]) for these diagrams yield 
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which is classified as improper or reducible in the nor- 
mal state£>^i2 and hence has been overlooked so far. 
This structure contains nothing that contradicts eq. (|5b[) , 
which defines self-energies, however. It is a natural con- 
sequence of the requirement that the Hugenholtz-Pines 
relation ([7]), which is an exact statement, is obeyed order 
by order in the self-consistent perturbation expansion. 
We will show shortly that it brings an additional term to 
eq. ([]}. On the other hand, the diagrams in the first and 
second rows in Fig. [5J which yield conventional proper or 
irreducible self-energies, are at least one order of magni- 
tude smaller at T — 0; thus, they can be neglected for 
the present purpose. 

The improper structure extends beyond the second or- 
der. In the third order, for example, the rightmost dia- 
gram in Fig. [3] also yields&ii 



£,(3ip) 
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where the factor 5/4 is the sum of 1 and 1/2 2 originating 
from the particle-hole and particle-particle bubble dia- 
grams, respectively. Equation (| 13[) reveals common fea- 
tures of the leading-order improper contribution at each 
order: (i) they are functionals of 



f-= G->+ F- — F- 



and (ii) the diagonal and off-diagonal self-energies are the 
same. Note the symmetry fp — fp. 

Adding the series of eq. (fl~3|) to eqs. ([8a)) and (|8bj) . we 
can express the diagonal and off-diagonal self-energies as 



Ap=A^ + A^\ (15) 



respectively, where A^ 
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satisfies 



A^ — A^ p \ The corresponding Green's functions are 



obtained by using eqs. (0 and (J7J as 
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This expression for fi appears to be the same as eq. (|8d[) , 
but the right-hand side now should be calculated by using 
eq. (|16j) instead of eq. ([9]) in the self-consistent approxi- 
mation scheme. To be specific, X^ 1 - 1 is still given by eq. 
U but with Gjp — > Gp in eq. 



cj, i.e., 



= «o -E4 ^ ,E,0+ - E( G ?- 4 )• ( 17b ) 



This redundant expression enables us to treat the third 
term on the right-hand side as a correction to the mean- 
field result (flbf . Similarly, eq. (|8b|) is replaced by 



A« = U 
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On the other hand, the substitution of eq. (fT6| into the 
series of eq. (fTB")) yields the equation for A^ = A~ lp ' + 

A (3ip) , 



as 



Aj p) = Uno 
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with 
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Equation (ITSl) forms an algebraic equation for A^ lp \ and 

ling ~ Un, i.e., it is of 



we see immediately that A^ 



(i P ) 



(14) the same order as A^ in eq. (j8b[) . 
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Let us transform the third term on the right-hand side 
of eq. (|17b|) by using eqs. © and p^|) , noting the sum- 
mation over p below eq. (f8c| at T = 0, and making a 
change of variables as 

€ p = Une p , si = Unii, Ap- = Un Ap. (20) 

The result may be written as 



where Cj p is given by 
4^ 



r 3/2 



depC^ 2 / cIq 
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(3 - 3)*?° 
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When substituting eq. (|2"Tj) into eq. (|17bp for the present 
purpose, we can replace (Un/8n) 3 ^ 2 by the leading-order 
expression (na) 3 / 2 with A^ — > 1 in eq. (|2"2"j) ; see eqs. 
dU, ©, and (Unj) on this point. Equation (jTHJ) is also 
approximated with Una ~ {/ ri as 



Ai 1 



(ip) _ 
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where 77^ is defined by 



(23) 



fj ff =(ej + e 2 p )/2e p +l. (24) 

Equation (23]) with eq. JM]) implies that Ai lp) is a func- 
tion of only x = £ 2 /e P + e P - Hence, it is convenient to 
make a further change of variables in eq. 



x=^ + e p , ^ = e p , b(x) = ^ + l, y(x) = A^\ 

(25) 

to transform it to 
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The last expression has been obtained by changing the 
order of integrations and subsequently integrating over 
£. A similar analysis on the third term in the square 
brackets of eq. (|17c[) yields 



(27) 



where cf is the same as eq. (|2"2")l except for e 2 — e 2 — > 
—e\ — e 2 in the numerator. Transforming it in the same 
manner as from eq. (|22j) into eq. (|26|) . we find that 



cf = 3c ip . 



(28) 



Let us substitute eqs. (J2TJ and (|2?J) with eq. J28J) into 
eq. (|17p . We then find that the improper contribution 
changes eqs. (fTbf and (|T2"j) into eqs. (|2b|) and 



/i = 8nan 



ip 
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respectively. Finally using the thermodynamic relation 
H = dE/dN to integrate eq. (|29l) over iV, we arrive at 
eq. (|2a| . Equation (|2|) with eq. (f26)) are correct up to the 
next-to-the-leading order. 

Finally, let us estimate c; p by using two different ap- 
proximations. First, retaining only the first two terms on 
the right-hand side of eq. (|23|) yields a cubic equation for 
y = y(x), which is defined in eq. (1251) . as 



y 



2by 2 + (b 2 + 2)y - 2(6 + 5) = 0. 



(30) 



It has a single real solution y = y(x) for < x < oo that 
approaches continuously as x — > oo, in accordance with 
A (ip )(ie<>) -s- as \e e \ -> oo. Substituting it into eq. (|26|) 
and performing the integration, we obtain c lp = 0.412. 
Hence, the factor 128/15-^ = 4.81 in eq. (|la[) is changed 
into (16/5)(8/3VtF + c ip ) = 6.13 in eq. (ga]). Second, 
eq. (|13b[) may be generalized to the nth-order improper 
terms (n = 3, 4, ■ • • ) originating from the particle-hole 
and particle-particle bubble diagrams as 
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where fp is defined by eq. (|T4"]) . Collecting the series for 
n = 3,4, ••• ,oo together with eq. (113a[) of the second 
order, we obtain 



I - 2Un fp 1 
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(26) This self-consistent equation for A^ is transformed in 
terms of the quantities in eq. ([23]) into 

y 4 - 3(6 - l)y 3 + (3b 2 - 66 + 4)y 2 
-(6 3 - 36 2 + 66 + A)y + 2{b 2 + 26 - 4) = 0, (33) 

which also has a single real solution y = y(x) for < x < 
oo that approaches continuously as x — > oo. Substitut- 
ing it into eq. (|26|) . we obtain q p = 0.563. Thus, we may 
expect that c; p is of order 1. The exact value of Ci P can 
only be reached by collecting all the leading-order im- 
proper terms, however, meaning that further studies are 
required for a quantitative estimate of c; p . On the other 
hand, a diffusion Monte Carlo study was performed on 
eq. ([I]); 13 the data in Table I of ref. for o?n < 10" 2 
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yield the value 4.3(1), rather than 128/15^ = 4.81 
from Eq. (|la|) . thus showing a clear deviation from the 
Lee-Huang- Yang expression beyond the numerical uncer- 
tainty. On the other hand, an experiment on eq. (jlal) in 
a trap potential reported a value 4.5(7)^ which still re- 
quires improvement for the present purpose, however. 
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